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Abstract 



This paper considers the class of stochastic processes X defined on [0, T] by X (t) — G {t, s) dM (s) 
where M is a square-integrable martingale and G is a deterministic kernel. When M is Brownian 
motion, X is Gaussian, and the class includes fractional Brownian motion and other Gaussian 
processes with or without homogeneous increments. Let m be an odd integer. Under the 
asumption that the quadratic variation [AI] of M is differentiable with E [\d [M] /dt\"^] finite, it 
is shown that the mth power variation 



In the case of a Gaussian process with homogeneous increments, 6 is X's canonical metric, 
the condition on S is proved to be necessary, and the zero variation result is extended to non- 
integer symmetric powers, i.e. using \X {s + e) — X (s)]™ sgn {X {s + e) — X (s)) for any real 
value m > 1. In the non-homogeneous Gaussian case, when m = 3, the symmetric (generalized 
Stratonovich) integral is defined, proved to exist, and its Ito formula is proved to hold for all 
functions of class C^. 

KEY WORDS AND PHRASES: Power variation, martingale Volterra convolution, co- 
variation, calculus via regularization, Gaussian processes, generalized Stratonovich integral, non- 
Gaussian processes. 
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1 Introduction 

The purpose of this article is to study wide classes of processes with zero cubic variation, and more 
generally, zero variation of any order. Before summarizing our results, we give a brief historical 
description of the topic of p- variations, as a basis for our motivations. 
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exists and is zero when a quantity 



(r) related to the variance of an increment of Af over a 




1 



1.1 Historical motivations 



The p- variation of a function / : [0, T] — > R is the supremum over all the possible partitions 
{0 = to < • • • < ^Af} of [0, T] of the quantity 

N-l 
i=l 

The analytic monograph [4] contains an interesting study on this concept, showing that a p- variation 
function is the composition of an increasing function and a Holder function. The notion of p- 
variation of a stochastic process or of a function was rediscovered in stochastic analysis, particularly 
in a context of pathwise (or quasi-pathwise) stochastic calculus. The fundamental paper [S], due 
to H. Follmer, treats the case of 2-variations. More recent dealings with p-variations and their 
stochastic applications, particularly to rough path and other integration techniques, are described 
at length for instance in the books [5j and [H] , which also contain excellent bibliographies on the 
subject. 

The Ito stochastic calculus has mimicked the notion of 2-variation, with the notion of quadratic 
variation. Given a semimartingale S, its quadratic variation [S] is the limit in probability of 

N-l 

Y,\S{ti+i)-S(^U)\^. (2) 
1=1 

For p = 2, instead of considering the supremum of quantities ([1]), one considers the limit of the 
expression ([2|) when the mesh of the partition converges to zero. Moreover, the notion becomes 
stochastic. In fact for a standard Brownian motion B, its 2-variation [B] is a.s. infinite, but its 
quadratic variation is equal to T. In order to reconcile 2-variations with the finiteness of [B], 
many authors have proposed restricting the supremum in ([T]) to the dyadic partitions. However, 
in Ito calculus, the idea of quadratic variation is associated with the notion of covariation (a.k.a. 
joint quadratic variation) [S^, 5^] of two semimartingales 5^, S"^, something which is not present is 
analytic treatments of p- variation. This covariation [S^, S'^] is obtained by polarization of ([2]), i.e. is 
the limit in probability of ^Z=~i^ {S^U+i) - S^U)) (S^iU+i) - S^{ti)) when, again, the partition 
mesh goes to zero. 

In the study of stochastic processes, the p-variation has been analyzed in some specific cases, 
such as local time processes (see [20j), iterated Brownian motion, whose 4-th variation is finite, and 
more recently fractional Brownian motion (fBm) and related processes. To work with a general 
class of processes, the tools of Ito calculus would nonetheless restrict the study of covariation 
to semimartingales. In [22J, the authors enlarged the notion of covariation to general processes 
X and Y. They achieved this by modifying the definition, considering regularizations instead of 
discretizations. An objective was to produce a more efficient stochastic calculus tool, able to go 
beyond the case of semimartingales. Given two processes X and Y, their covariation [X, Y] is the 
limit in probability, when e goes to zero, of 

(X{s + e)-X{s))(Y{s + e)-Y{s)) 

C{e,X,Y){T)= ^ "lJI^ i L^ds; t>0. (3) 

Jo £ 

The limit is again denoted by [X, Y] and this notion coincides with the classical covariation when 
X, Y are continuous semimartingales. The processes X such that [X, X] exists are called finite 
quadratic variation processes; their analysis and their applications were performed in [7} I21j. 
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The notion of covariation was also extended in order to involve more processes. In [lOj the 
authors consider the n-covariation of n processes X^,... ,X^. For n = 4, for X being an fBm 
with so-called Hurst parameter H = 1/4, the paper [S] calculates the covariation [g (X) ,X,X,X] 
where g is, for instance, a bounded continuous function. If X = X^ = X'^ = X^ is a single 
stochastic process, we denote [X;3] := [X, X, X], which is called the cubic variation, and is one of 
the main topics of investigation in our article. Note that this variation involves the signed cubes 
{X (s + e) — X(s))^, which has the same sign as the increment X {s + e) — X(s), unlike the case of 
quadratic or 2-variation, or of the so-called strong cubic variation, where absolute values are used 
inside the cube function. The significance of the cubic variation lies in its ability to guarantee the 
existence of (generalized symmetric) Stratonovich integrals, and their associate Ito formula with no 
correction term, for highly irregular processes, notably fBm with if > 1/6 (such fBm have Holder 
regularity parameter exceeding 1/6; compare with the near 1/2-Holder-regularity for continuous 
semimartingales) . This is a relatively well-known phenomenon, established in great generality in 
the landmark paper [12]. This result is a main motivation for our work, in which we attempt to give 
specific classes of Gaussian and non-Gaussian processes with the said zero cubic variation property; 
indeed the conditions in pL2j have really only been applied to the very specific Gaussian class of 
fBm, which are fractionally self-similar and have stationary increments. 

1.2 Summary of results 

Let [0,T] be a fixed time interval. Following the regularization methodology of [22j or [23j . the 
cubic variation of a process X, denoted by [X, 3], was defined similarly to [TO] as the limit in 
probability or in the mean square, as e ^ 0, of {X{s + e)-X {s)f ds. This wih be nuh for 

a deterministic function X as long as it is a-Holder-continuous with a > 1/3. But the main physical 
reason for being interested in this cubic variation for random processes is that, because the cube 
function is symmetric, if the process X itself has some probabilistic symmetry as well (such as the 
Gaussian property and the stationarity of increments), then we can expect [X,3] to be for much 
more irregular processes than those which are almost-surely a-Holder-continuous with a > 1/3. As 
mentioned above, it is known from [TO] that fBm, which is parametrized by the Hurst index H, for 
which one has that X is a-Holder-continuous almost surely for all a > H, has zero cubic variation for 
all H > 1/6. This doubling improvement over the deterministic situation is due exclusively to the 
random symmetries of fBm, as they combine with the fact that the cube function is odd. Typically 
for other types of variations, such as the strong cubic variation \ X (s + e) — X {s)\^ ds, H 

needs to be larger to guarantee existence of the variation, let alone nullity. 

These properties of fBm beg the question of what occurs for other Gaussian processes which 
may not be self-similar or even have stationary increments, or even for non-Gaussian processes with 
similar Holder-continuous paths, and to what extent the threshold a > 1/6 is sharp. Similarly, 
can the odd symmetry of the cube function be generalized to any "symmetric" power function, i.e. 
X i—s- |x|™'sgn(x) with arbitrary integer or non-integer m > 1 ? This refers to what we will call the 
"odd mth variation" , defined (when it exists in the mean-square sense) by 

[X, m] := lim e'^ [ \X {s + e) - X (s)r sgn {X {s + e) - X (s)) ds. 
Jo 

The qualifier "odd" above, when applied to m = 3, can easily yield the term "odd cubic variation", 
which has historically been called simply "cubic variation" as we noted before, as opposed to the 
"strong cubic variation" which involves absolute values; therefore in this article, we will systemat- 
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ically use the qualifier "odd" for all higher order mth variations based on odd functions, but will 
typically omit it for the cubic variation. 

This article provides answers to some of the above questions, both in Gaussian and non-Gaussian 
settings, and we hope that it will stimulate work on resolving some of the remaining open problems. 
Specifically, we consider the process X defined on [0, T] by 



X{t)= [ G {t, s) dM (s) (4) 
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where G is a non-random measurable function on [0, T] , which is square-integrable with respect 
to s G [0,7"] for every fixed t, and M is a square-integrable martingale on [0, T]. In other words, 
X is defined using a Volterra representation with respect to a square-integrable martingale. The 
quadratic variations of these martingale-based convolutions was studied in [6]. 

What we call the Gaussian case is that where M is the standard Wiener process (Brownian 
motion) W. The itemized list which follows is a summary of our results. There are several technical 
conditions which are needed in our various result statements, and are increasingly restrictive as we 
move from simple Gaussian cases to non-Gaussian cases, but nonetheless represent wide classes 
of processes. Here, for the reader's convenience, we have not spelled out all these conditions, 
indicating instead references to the precise theorem statements in the body of this article. One 
condition which appears in all cases, is essentially equivalent to requiring that all processes X 
that we consider are not more regular than standard Brownian motion, i.e. are not 1/2-Holder- 
continuous. This typically takes the form of a concavity condition on the process's squared canonical 
metric 5^ (s, t) := E (X (t) — X (s))^ . This condition is hardly a restriction since the main interest 

of our results occurs around the Holder exponent 1/6, or more generally the exponent l/(2m) for 
any m > 1: the processes with zero odd mth variation appear as those which are better than 
1 / (2m)-H61der-continuous. 

• [Theorem [6] on page II Ij. When X is Gaussian with homogeneous increments (meaning 6{s,t) 
depends only on |t — s|), for any odd integer m > 3, X has zero odd mth variation if and only 
if 5(r) = o(ri/(2'")). 

• [Theorem [9] on page [19]. The sufficient condition of the result above holds for any real m > 1. 

• [Theorem [7] on page [16]. When X is Gaussian with non-homogeneous increments, for any 
odd integer m > 3, if 6"^ {s,s + r) = o (r^/^^™)) uniformly in s, and under a technical non- 
explosion condition on the mixed partial derivative of 5^ near the diagonal, X has zero odd 
mth variation. 

• [Theorem [To] on page I23j. Let m > 3 be an odd integer. When X is non-Gaussian as in (j4]), 
based on a martingale M whose quadratic variation process has a derivative with 2m moments 
(the actual condition on M in the theorem is weaker), let F (t) = (E [{d [M] /dt)"'])^/^^"') and 
consider the Gaussian process 



Z{t)= [ r{s)G{t,s)dW (s) 
Jo 



Under a technical integrability condition on planar increment of TG near the diagonal, if Z 
satisfies the conditions of Theorem [6] or Theorem [71 then X has zero odd mth variation. 
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Proposition [TT] on page [27] provides the example of a wide class of martingales and kernels for 
which the assumptions of Theorem [10] are satisfied. 

[Theorem [12] and its corollary, on pagel31|. When m > 3 is an odd integer and X is a Gaussian 



process with non-homogeneous increments such that 5 {s,s + r) = o 



a/{2m)^ 



uniformly in 



s, under some additional technical conditions whose scope is discussed after the corollary on 
page [321 foi' every bounded measurable function g on R, 



lim — 

£^0 



E 



du {Xu^ 



0. 



[Corollary [T3l on page [32] . If m = 3, by results in [12j . the theorem implies that for any 
/ G (R) and t > 0, the Ito formula / {Xt) = f (Xq) + f X (odu) holds, where the 
integral is in the symmetric (generalized Stratonovich) sense. 



1.3 Relation with other recent work 

We finish this introduction with a description of recent work done by several other authors on 
problems related to our preoccupations to some extent, in various directions. The authors of the 
preprint [13] consider, as we do, stochastic processes which can be written as Volterra integrals with 
respect to martingales, just as we did. In fact, they study the concept of "fractional martingale", 
which is the generalization of the so-called Riemann-Liouville fractional Brownian motion when 
then driving noise is a martingale. This is a special case of the processes we consider in Section 
[U with K(t,s) = {t — s)^~^^^. The authors' motivation is to prove an analogue of the famous 
characterization of Brownian motion as the only continuous square-integrable martingale with a 
quadratic variation equal to t. They provide similar necessary and sufficient conditions based on 
the 1///- variation for a process to be fractional Brownian motion. The paper [13j does not follow, 
however, the same motivation as our work: for us, say in the case of m = 3, the threshold H > 1/6 
is the key point for (odd) cubic variations, with its relation to stochastic calculus. 

To find a similar motivation as ours, one may look instead at the equally new result of [17] . 
where the authors study the central and non-central behavior of weighted power variations for fBm. 
Using the Hermite polynomial of order q rather than the power-g function, largely for simplicity 
(we do not believe that this is a necessary restriction, and we think that the authors would agree 
with us), they show that the threshold value H = 1/ {2q) poses an interesting open problem, since 
above this threshold (and below H = 1 — l/{2q)) one obtains Gaussian limits (these limits are 
conditionally Gaussian when weights are present, and can be represented as stochastic integrals 
with respect to an independent Brownian motion), while below the threshold, degeneracy occurs, 
and the behavior at the threshold is not known except (as far as we understand the state of the 
art at the time of writing this) for g' = 2,3. In [17J, a threshold Ito formula for q = 3, i.e. at 
H = 1/6, boasts an intriguing exotic correction term with an independent Brownian motion, and 
raises important questions; for instance, one can compare this with a similar result in [16| for 
H = 1/4, but this time with bidimensional fBm, in which two independent Brownian motions are 
needed to characterize the exotic correction term. These questions extend beyond the case of fBm, 
to our framework, beyond the scope of our article. 

The value H = 1/6 is mentioned again in the context of the stochastic heat equation driven 
by space-time white-noise, in which discrete trapezoidal sums converge in distribution (not in 
probability) to a conditionally independent Brownian motion: see [3] and [18]. In fact, one can 
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show for fBm with H = 1/6, that the cubic variation in the regularized sense does not exist in 
the sense of convergence in probabihty. Specifically, using arguments developed in [11], the paper 
[T2] proves in Theorem 4.1 part (2), convergence in distribution to a normal. Using some recently 
established arguments about the convergence of sequences in fixed Wiener chaoses, one can then 
show (see Proposition [5] herein) that this convergence cannot be in probability. By establishing a 
necessary and sufficient condition around the same threshold for nullity of the cubic variation for 
many Gaussian processes with homogeneous increments, our paper confirms the value H = 1/6 as 
a crucial regularity parameter for existence of cubic variations, showing that self-similarity is not 
related to such existence. 

Summarizing, when compared to the works described in the above three paragraphs, our work 
situates itself by 

• choosing to prove necessary and sufficient for nullity of the cubic variation (around the thresh- 
old regularity value H = \/&) for Gaussian processes with homogeneous increments (a wider 
class than previously considered), and studying this threshold for higher "odd" power func- 
tions, with possibly non-integer order; 

• showing that our method is able to consider processes that are far from Gaussian and still 
yield sharp sufficient conditions for nullity of odd power variations, since our base noise may 
be a generic martingale with only a few moments. 



The article has the following structure. Section [2] contains some formal definitions and notation. 
The basic theorems in the Gaussian case are in Section [3l where the homogeneous case, non- 
homogeneous case, and case of non-integer m are separated in three subsections. The use of 
non-Gaussian martingales is treated in Section [H Section [5] presents the Ito formula. 



2 Definitions 

We recall our process X defined for all i G [0, by 

X{t)= ! G{t,s)dM{s) (5) 

^0 

where G is a non-random measurable function on [0, T]^, which is square- integrable with respect to 
s G [0, T] for every fixed t, and M is a square-integrable martingale on [0, T]. For any real number 
m > 2, let the odd e-m-th variation of X be defined by 

[X,m],{T):=- [ ds\X{s + e)-X{s)rsgn{X{s + e)-X{s)). (6) 
£ Jo 

The odd variation is different from the absolute (or strong) variation because of the presence of the 
sign function, making the function |x|™'sgn(x) an odd function. In the sequel, in order to lighten 
the notation, we will write (x)"^ for |x|'"sgn(x). We say that X has zero odd m-th variation (in 
the mean-squared sense) if the limit 

lim[X,m]aT) = (7) 

holds in (0). 
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The canonical metric 6 of a stochastic process X is defined as the pseudo-metric on [0, T]^ given 



by 



d^{s,t) = 'E {Xit)-Xis)y 



The covariance function of X is defined by 

Q{s,t)=^[X {t)X{s)]. 

The special case of a centered Gaussian process is of primary importance; then the process's entire 
distribution is characterized by Q, or alternately by 5 and the variances var {X (t)) = Q (t, t), since 
we have Q{s,t) = ^ {Q {s, s) + Q {t,t) — 6 {s,t)). We say that 6 has homogeneous increments if 
there exists a function on [0, T] which we also denote by 6 such that 

S{s,t) = 6i\t-s\). 

Below, we will refer to this situation as the homogeneous case. This is in contrast to usual usage of 
this appellation, which is stronger, since for example in the Gaussian case, it refers to the fact that 
Q {s, t) depends only on the difference s — t; this would not apply to, say, standard or fractional 
Brownian motion, while our definition does. In non-Gaussian settings, the usual way to interpret 
the "homogeneous" property is to require that the processes X [t + ■) and X (•) have the same law, 
which is typically much more restrictive than our definition. 

The goal of the next two sections is to define various general conditions under which a charac- 
terization of the limit in ([7]) being zero can be established. In particular, we aim to show that X 
has zero odd m-th variation for well-behaved M's and G's as soon as 



5{s,t) = o{\t 



l/(2m)^ ^ 



(8) 



and that this is a necessary condition in some cases. Although this is a mean-square condition, it can 
be interpreted as a regularity (local) condition on X] for example, when X is a Gaussian process with 
homogeneous increments, this condition means precisely that almost surely, the uniform modulus of 
continuity w of X on any fixed closed interval, defined by w (r) = sup{|X (t) — X (s)| ■.\t — s\ < r}, 

satisfies cj (r) = o ^r^/^ log^^^ (I/'")) • 

3 Gaussian case 

We assume that X is centered Gaussian. Then we can write X as in formula ([5]) with M = W a 
standard Brownian motion. More importantly, beginning with the easiest case where m is an odd 
integer, we can easily show the following. 

Lemma 1 If m is an odd integer > 3, we have 



E 



{[X,m]eiT)f 

(m-l)/2 T i-t 

^ / dtds& {s,tr-^^ 

^ f^o Jo Jo 

(m-l)/2 
j=0 



Var [X{t + e)-X {t)y Var [X {s + e) - X (s) 
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where the cj 's are constants depending only on j, and 

{s, t) := E [{X {t + e)-X (t)) (X {s + e) - X (s))] . 

Proof. The lemma is an easy consequence of the fohowing formula, which can be found as 
Lemma 5.2 in |12j : for any centered jointly Gaussian pair of r.v.'s {Y, Z), we have 

(m-l)/2 

E [y^Z'"] = CjE [YZ]""-^^ Var [X^ Var [Y]> . 

j=0 

m 

We may translate is,t) immediately in terms of Q, and then 5. We have: 

{s,t) = Q {t + e,s + e) - Q (t, s + e) - Q {s,t + e) + Q {s, t) 

= \ [-^'^ {t + £,s + e) + 5'^ {t, s + e) + 6^ {s,t + e) - 6^ (s, t)] (9) 

= : -2^is,tyXs+e,t+e)^'^- (10) 

Thus {s,t) appears as the opposite of the planar increment of the canonical metric over the 
rectangle defined by its corners (s, t) and {s + e,t + e). 

3.1 The case of fBm 

Before finding sufficient and possibly necessary conditions for various Gaussian processes to have 
zero cubic (or mth) variation, we discuss the threshold case for the cubic variation of fBm. Recall 
that when X is fBm with parameter H = 1/6, as mentioned in the Introduction, it is known from 
[12\ Theorem 4.1 part (2)] that [X, 3]^ (T) converges in distribution to a non-degenerate normal law. 
However, there does not seem to be any place in the literature specifying whether the convergence 
may be any stronger than in distribution. We address this issue here. 

Proposition 2 Let X be an fBm with Hurst parameter H = 1/6. Then X does not have a cubic 
variation (in the mean-square sense), by which we mean that [X, 3]^ (T) has no limit in Lp' (i7) as 
e ^ 0. In fact more is true: [X, 3]^ (T) has no limit in probability as e — > 0. 

In order to prove the proposition, we study the Wiener chaos representation and moments of 
[X, 3]e (T) when X is fBm; X is given by ([5]) where W is Brownian motion and the kernel G is 
well-known. Information on G and on the Wiener chaos generated by W can be found respectively 
in Chapters 5 and 1 of the textbook [19]. The covariance formula for an fBm X is 

Rh {s, t) := E [X (t) X (s)] = 2'^ (s^// ^ ^2H _ |^ _ ^^2H^ _ ^^^^ 
Lemma 3 Fix e > 0. Let AXs := X {s + e) - X {s) and AGs (u) := G {s + e,u) - G {s,u). Then 



[X,3],(T) =Ti+T; 



3 



ds / AGsiu)dW{u) / \AGs{v)\'^dv (12) 



3 

£ Jo Jo ' 'Jo 

6 



_ dwiss) r dw{s2) r dw{si) r 

£ Jo Jo Jo Jo 



n {sk) 

.k=l 



ds. (13) 
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Proof. The proof of this lemma is elementary. It follows from two uses of the multiplication 
formula for Wiener integrals [19\ Proposition 1.1.3], for instance. It can also be obtained directly 
from Lemma 13.21 below, or using the Ito formula technique employed further below in finding an 
expression for [X, mje (T) in Step of the proof of Theorem [TOl on page[23l All details are left to 
the reader. ■ 

The above lemma indicates the Wiener chaos decomposition of [X, 3]^ (T) into the term Ii of 
line (jl2p which is in the first Wiener chaos (i.e. a Gaussian term), and the term T3 of line (jl3p . in 
the third Wiener chaos. The next two lemmas contain information on the behavior of each of these 
two terms, as needed to prove Proposition [2l 

Lemma 4 The Gaussian term Zi converges to in {n) as e ^ 0. 

Lemma 5 The 3rd chaos term T3 is bounded in ($7) for all e > 0, and does not converge in 
L2 (n) as e^O. 

Proof of Proposition [2], We prove the proposition by contrapositive. Assume [X, 3]^ (T) 
converges in probability. For any p > 2, there exists Cp depending only on p such that E [l^il^] < 



E 



)''^'andE[|J3n <Cp(E jXal 



p/2 



; this is a general fact about random variables in 



fixed Wiener chaoses, and can be proved directly using Lemma [3] and the Burkholder-Davis-Gundy 

+E 



) < 00 by Lemmas m and [5l we also 
ay uniform integrability, [X, 3]^ (T) = Ti +X3 



inequalities. Therefore, since we have suP(,>q(E 
get sup£^q(E < 00 for any p. Therefore, 

converges in (Q). In (J7), the terms Xi and I3 are orthogonal. Therefore, Ii and T3 much 
converge in (Q,) separately. This contradicts the non-convergence of X3 in (Q,) obtained in 
Lemma [5l Thus [X, 3]^ (T) does not converge in probability. ■ 

To conclude this section, we only need to prove the above two lemmas. To improve readability, 
we write H instead of 1/6. 

Proof of Lemma |4l From Lemma [1] with m = 3, and Lemma [3l we get 



E 



-e 
2 



2' 



JQ 



'2+4H 



-2+4H 



dtdsQ' {s, t) Var {X {t + e) - X (t)) Var {X {s + e) - X (s)) 
/ / dtds Cov [X {t + e) - X it) -X {s + e) - X (s)] 

JQ Jo 

[ {X{t + e)-Xit)), 
Jo 



Var 



Var 



T+e 



X (t) dt 



\dt 



X (t) dt 



Bounding the variance of the difference by twice the sum of the variances, and using the fBm 
covariance formula (1111). 



E 



< 9e 



< 9e 



-2+4H 



'2+AH 



T+e rT+e 



Rh {s, t) dsdt + Rh (s, t) dsdt 

Jo Jo 



proving Lemma [H 
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Proof of Lemma [5l From Lemma [T] with m = 3, Lemma [3l and the covariance formula (jlip 
for fBm, we get 



E 



T rt 



12 

£ Jo Jo 

1' 

Jo 



dtds {Q^{s,t)f 



dtds [\t-s + e\^" + \t-s- e\^" -2\t- sp^ 



We must take care of the absolute values, i.e. of whether e is greater or less than t — s. We define 
the "off-diagonal" portion of E 



3 •= 



6e" 



2e Jo 



as 



dtds {\t-s + ep^ + \t-s- - 2 It - 



For s,t in the integration domain for the above integral, since t := t — s > 2e, hy two iterated 
applications of the Mean Value Theorem for the function x"^^ on the intervals [i — e, i] and [t, t + e], 



|t + ep^ + \t-e\ 



2H _ ^^H ^ 2H {2H - 1) e (a - 6) e 



2H-2 



for some ^2 S [t — e, t], ^1 G [t, t + e], and ^ G [Cij ^2], and therefore 



lOPJsl < 24H\2H - l|e~2 



T rt~2e 



2AH\2H - We'^ 

l2e Jo 

2AH\2H - 11 /"^ r 



2e JO 
T ft-2e 



\2H-2 



5-6H 



2e 



e{t — s — e 
{t-s- ef"-" dtds 



dtds 



2m\2H-l\ fifT 1 2AH\2H-l\ 

< ^ [j, m~l ^ ^ [rp 

- 5-6H 5-6H 
Thus the "off-diagonal term is bounded. The diagonal part of X3 is 



3 • = 



6e" 



Jt~2e 
2e 



dtds {\t-s + e\^^ + \t-s- e|^^ -2\t- sp^ 



di ( It + e|2^ + It - el^-^ - 2 ItT^ 



6e 



-l+6Hr 



dr (\r 



+ l|^-^ + |r-l|^'' -2|r|^''j dr = C{H)T 



\2H 



2H\ 



where C (H) only depends on H. Thus the diagonal part PT3 of E[|23p] is constant. This proves 
that X3 is bounded in (O), as announced. To conclude that it cannot converge in ($7), recall 
that from [12, Theorem 4.1 part (2)], [X, Sj^ (T) = Xi + X3 converges in distribution to a non- 
degenerate normal law. By Lemma HI Xi converges to in {n). Therefore, X3 converges m 
distribution to a non-degenerate normal law; if it also converged in (Q), since the 3rd Wiener 
chaos is closed in (17), the limit would have to be in that same chaos, and thus would not have 
a non-degenerate normal law. This concludes the proof of Lemma [5l ■ 
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3.2 The homogeneous case 

We now specialize to the homogeneous case. We are ready to prove the first sufficient condition for 
having a zero m-th variation when m is an odd integer. 

Theorem 6 Let m > 1 be an odd integer. Let X be a centered Gaussian process on [0, T] with 
homogeneous increments; its canonical metric is 



5^is,t) :=E iX{t)-X{s)f 



5^(\t-s\ 



where the univariate function d'^ is increasing and concave on [0, T]. Then X has zero mth variation 
if and only if 5 (r) = o (^r^/C^'^)^ 

Proof. Step 0: setup. We denote by d6'^ the derivative, in the sense of measures, of 5^; we 
know that d6'^ is a positive bounded measure on [0, T]. Using homogeneity, we also get 

Var [X{t + e)-X (t)] = 6^ (e) . 

Using the notation in Lemma [H we get 

rT ft 

\m-2j 



jj = £-2 ^, f dt [ dse' (s, ty 

Jo Jo 



Step 1: diagonal. Let us deal first with the diagonal term. We define the e-diagonal :- 
{0 < t — e < s < t < T}. Trivially using Schwartz's inequality, we have 



is,t)\ < ^Yar {X{t^e)-X (t)\ Var [X {s + e) - X (s)] = 5^ (e) 



Hence, according to Lemma [U the diagonal portion ^j^Q Jj,D of E {[X,m]^ {T))^ 
bounded above, in absolute value, as: 



can be 



(m-l)/2 

E 



(m-l)/2 



E '^^'^■■= E ^-'^'ne)cJ dt! 

^ (m-l)/2 



dsG^ {s,t) 



m-2j 



j=0 



dsS''^ (e) = est ■ e-i<52™ (e 



The hypothesis on 5'^ implies that this converges to as e tends to 0. 

Step 2: small t term . The term for t £ [0, e] and any s £ [0, t] can be dealt with similarly, and is 
of a smaller order than the one in Step 1. Specifically we have 



\J- 



j,s\ ■■= e 



e Co 



dt / dse^s,t) 
Jo 



m—2j 



< (e) Cj(52(™-2J) (e) = ^.^^m ^ 



which converges to like o(e). 
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Step 3: off-diagonal. Because of the homogeneity hypothesis, we can calculate from Q that for 
any s, t in the e-off diagonal set OD^ := {0 < s < t — e < t < T} 



{s, t) = [5'^{t-s + e)- 52 {t - s)) - [5^ it-s)-6^{t-s- e)) 



t—s+e 



(16^ (r) 



t-s 



dS"^ (r) 



(14) 



t—s—e 



By the concavity hypothesis, we see that 0^ (s,t) is negative in this off-diagonal set OD^. Unfor- 
tunately, using the notation in Lemma [H this negativity does not help us because the off-diagonal 
portion Jj^oD of Jj also involves the constant cj, which could itself be negative. Hence we need to 
estimate Jj^oD more precisely. 

The constancy of the sign of 0^ is still useful, because it enables our first operation in this step, 
which is to reduce the estimation of {Jj^odI to the case of j = (m — 1) /2. Indeed, using Schwartz's 
inequality and the fact that |0^| = —0^, we write 



\Jj,OD\ = £ 



< {e) 



e Co 



dt 

T 

dt 



t-e 



ds\e' {s,t)\ 



m—2j 



dt 



t-e 



dsQ' (s,i)|0" (s,t)| 



m-2j-l 



ds{-Q' {s,t)) (e) 



m— 2j — 1 



= £-^5^"^-^ {£)\cj\ J dt J ds{-Q'{s,t)). 
It is now sufficient to show that the estimate for the case j = (m — 1) /2 holds, i.e. that 



dt 



t-e 



ds{-e' {s,t)) < ceb^ (2e) 



(15) 



We rewrite the planar increments of (5^ as in (jl4p to show what cancellations occur: with the 
notation s' = t — s, 

-e' (s, t) = - {S^ {s' + e)-6^ {s')) + {6^ {s') - 6^s' - e)) = - f d5^ (r) + d5^ (r) . 

Js' Js'-e 

Therefore, using the change of variables from s to s', and another to change [s' — £, s'] to [s', s' + e], 



dt 



t-e 



ds{-Q' {s,t)) 



dt 



dt 



dt 



dt 



ds' / dJ^ _ / 



ds' / dS"^ (r) - / ds' d5^ (r) 



t-e 



ds" / d5^ (r) - / ds' d5^ (r) 

Js" Je Js' 

rs"+e f-t rs'+e 

ds" / dd"^ (r) - / ds' d5^ (r) 

Js" Jt-e Js' 



(16) 
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We may now invoke the positivity of d6'^, to obtain 



dt 



ds{-Q^s,t)) < dt ds" 



d6^ (r) 



= / dt r ds" {5^ {s" + e)-5^ {s")) < [ die 6^ (2e) < Te6^ {2e) . 

Je Jo Je 

This is precisely the claim in (jlSp . which finishes the proof that for all j, \ Jj^oD\ ^ ce~^6'^"^ (2e) for 
some constant c. Combining this with the results of Steps 1 and 2, we obtain that 



E 



i[X,m]eiT)y <cst-£-'6'"^{2e) 



which implies the sufficient condition in the theorem. 

Step 4^ necessary condition. The proof of this part is more delicate than the above: it requires an 
excellent control of the off-diagonal term, since it is negative and turns out to be of the same order 
of magnitude as the diagonal term. We spell out the proof here for m = 3. The general case is 
similar, and is left to the reader. 

Step positive representation. The following lemma is proved in the appendix. It uses the 
following chaos integral notation: for any n G N, for g ^ ([0, T]"), g symmetric in its n variables, 
then In (g) is the multiple Wiener integral of g over [0, T]" with respect to W. 

Lemma Let f e ([0, T]). Then h iff = 3 |/li2([o,T]) h if) + h{f^f^f) 

Using this lemma, as well as definitions ([5]) and ([6|), and introducing the notation AG, (u) = 
G {s + e,u) — G (s, u), and exploiting the fact that the covariance of two multiple Wiener integrals 
of different orders is 0, we can write 



E 



4 / ds 



1 

£ Jo 



T 



-J2 i 



dtE 



{X{s + e)-X {.s)f iX{t + £)-X {t)Y 



dtE 



hiAGsfhiAGtf 
ds dtE [h [AGs) h (AGt)] I AG,|i2([o_y]) | AGtl^^ 



([0,T]) 



+ ds 



T 



dtE 



h({AGsfnh[{AGt) 



Now we use the fact that E [I3 (f) I3 (g)] = {f ^9) l'^{\qt]^)i plus the fact that in our homogeneous 



situation |AGs|^2 



{[0,T]) 



([0,T]3). 

52 (e) for any s. Hence the above equals 



L2([0,T]3) 



ds dt{AGs,AGt)LH\oT]) + - ds ( (AG,)®^ (AG^ 

£ Jo Jo U : J/ £ \ 

duAGs (n) AGt (u) + ^ ds '^^JJJ'^ {du^AG^ {ui) AGt (n^)) 



(e) 



ds 



dt 



[0,T]3 



i=l 



%^ ie) 



du 



dsAGs (u) 



^ J J '^'"1 dU3 
[0,T]3 




dslliAG, M) 
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Step 4-2: Ji as a lower bound. The above representation is extremely useful because it turns out, 
as one readily checks, that of the two summands in the last expression above, the first is what we 
called Ji and the second is Jq, and we can now see that both these terms are positive, which was 
not at all obvious before, since, as we recall, the off-diagonal contribution to either term is negative 
by our concavity assumption. Nevertheless, we may now have a lower bound on the e-variation by 
finding a lower bound for the term Ji alone. 

Reverting to our method of separating diagonal and off diagonal terms, and recalling by Step 2 
that we can restrict t > e, we have 



qxi f-T ft fT 

Ji = ^2 / dt ds duAGs (u) AGt (u) 

e Je Jo Jo 



9<5* (e) 



2 [ dt I dsQe{s,t) 

£ Je Jo 



nxA (p\ f-T ft 

^ ' ' dt ds{6^{t-s + e)-6^{t-s)- {6^t-s)-5^{\t-s-e\))) 



£ Je JO 
= Jl,D + Jl,OD 

where, performing the change of variables t — s i— > s 

J, ,, := r dt f ds {6^ {s + e)- 6^ (s) - {6^ (s) - 6\e - s))) 



£ Je Jo 

j^od:=^-^^I dt I ds[5\s + e)-6\s)-{5\s)-6^{s-e))). 



Step 4-3: Upper bound on |Ji,od|- Using the calculations performed in Step 3 (note here that 
(m — 1) /2 = 1, in particular line ()16p . we have 



Ji,OD = — I dt 



QA4 (c-^ rT r /■* rs+e re /-s+e 



e 



ds I d6' (r) - I ds I d6^ (r) 

t-e 



--■.K1+K2. 



We can already see that -fCi > and K2 < 0, so it's only necessary to find an upper bound on 
\K2\; but in reality, the reader will easily check that \Ki\ is of the order 5^ (e), and we will see that 
this is much smaller than either Ji^o or \K2\. Performing a Fubini on the variables s and r, the 
integrand in K2 is calculated as 



fe fs+e fe fs+e fe fr fZe fe 

/ / (^) = / / ^ / ^^2 / ^^2 / 

Jo Js Jo Js Jr=0 Js=0 Jr=e J s=r-e 

re r2e 

r dS"^ (r) + / (2e - r) dd"^ (r) 



r=0 



[r52 (^)] J _ [^^2 (^)j _ ^2 ^2 ^ 2e (5^ (2^) - 5^ (e)) 





e f2e 
5^ (r) dr + I 5^ (r) dr. 
Je 



2e 
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In particular, because \Ki \ <^ \K2\ and S is increasing, we get 

^4 pT / /•2£ 



\Ji,Od\ < ^ dtQ'^^ 6' (r) dr - 6' (r) dr) 



S'fep Lower bound on J\^d- Note first that 

/ ds ((^2 (s) - ^2 _ g-)^ = [ ds5^ (s) - / (5^ (e - s) = 0. 
7o Jo Jo 



Therefore 



9<5'i (e) 



[ dt [ ds ((5^ (s + £) - (5^ (s)) 

9^4 re rs+e 



£2 



(T-e) f ds f dd^ (r) 



We can also perform a Fubini on the integral in Ji^d, obtaining 



/•£ rs+e re rze 

ds dS^ {r)= r dS^ (r) + e / dS^ (r) 

Jo J s Jo Je 

= [rS^ (r)] I- (^) dr + e {5^ (2£) - 6^ (e)) 

= £(5^ {2e) - [ (r) dr. 
Jo 



In other words, 



£■ 



Ji,i. = ^^(r-£)(£(5M2£)- / 5^(r)cir). 



Step 4-5: conclusion. We may now compare Ji^d and |Ji,oDh using the results of Steps 4.1 and 
4.2, 



-/i — — Ki,o£>l ^ — — 



(T - £) (^£(52 (2£) - 6^ (r) dr j 



£2 



(T ^'^ ir) dr - 5'^ (r) dr^ 

(T - £) (2£) _ 52 (r)) dr. 



£2 



When 5 is in the Holder scale 5 (r) = r^, the above quantity is obviously commensurate with 
8^ {e) /£, which implies the desired result, but in order to be sure we are treating all cases, we now 
present a general proof which only relies on the fact that 5"^ is increasing and concave. 
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Below we use the notation ((5^)' for the density of which exists a.e. since 5^ is concave. 
Throwing away the integral from e to 3e/2, using the mean value theorem and the monotonicity of 
(5^ and ((5^) , we can write 

Ji > est ■ (e) (^5^ (2e) - 6^ (^2e^ 

> est -8-^6^ (e)-(e/2)- (52)'(2e). 

Since 5'^ is concave, and 5 (0) = 0, we have 6'^ (e) > 5^ (2e) /2. Hence, with the notation / (x) = 
6^ (2x) , we have 

Ji>cst-f (e) f {e) = cst-{fy{e). 

Therefore we have that lim^^o {f) (e) = 0. We prove this implies lim£_>o£ ^/^ (e) = 0. Indeed, 
fix r/ > 0; then there exists > such that for all e G (0,e,,], < {f^)' (e) < '>] (we used the 
positivity of {S"^)')- Hence, also using / (0) = 0, for any e G (0, e^], 

0<Zl(£l = i r (f)' {x)dx <- rr]dx = rj. 
^ ^ Jo ^ Jo 

This proves that lim^^o^^^/^ (s) = 0) which is equivalent to the announced necessary condition, 
and finishes the proof of the theorem. ■ 

3.3 Non-homogeneous case 

The concavity and homogeneity assumptions were used heavily above for the proof of the necessary 
condition in Theorem [6l However, these assumptions can be considerably weakened while still 
resulting in a sufficient condition. We now show that a weak uniformity condition on the variances, 
coupled with a natural bound on the second-derivative measure of 5^, result in zero m- variation 
processes. 

Theorem 7 Let m > 1 be an odd integer. Let X be a centered Gaussian proeess on [0, T] with 
canonical metric 

6^ {s,t) :=E [(X {t)-X{s)y 
Define a univariate function on [0, T], also denoted by 5"^, via 

6"^ (r) := sup 5^ (s, s + r) , 

se[o,T] 

and assume that 

5(r) =o(ri/2m^ . (18) 

Assume that, in the sense of distributions, the derivative dd"^ / (dsdt) is a finite signed measure /_i 
on [0, T]^; denote the off-diagonal simplex by OD = {{s,t) ■.0<s<t — e< T}; assume ^ satisfies, 
for some constant c and for all e small enough, 

1^1 (OD) < ce-^+^/™, (19) 

where is the total variation measure of ji. Then X has zero mth variation. 
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Proof. Step 0: setup. Recall that by Lemma [H 

(m-l)/2 J, . 



E 



i[X,m]eiT)y 



f2 J2 

(m-l)/2 
j=0 



JO 



dtdsQ' {s,tr-^^ 6^^{s,s + e) 6^^t,t + e). (20) 



and now we express 



9''(s,t) =^([s,s + e] X [t,t + e)) = / / ^{dudv) 

Js Jt 



(21) 



We again separate the diagonal term from the off-diagonal term, although this time the diagonal 
is twice as wide: it is defined as {(s, t) : < t — 2e < s < t}. 

Step 1: diagonal. Using Schwartz's inequality which implies |0^ (s, t)\ < 5 {s ^ s + e) 5 {t,t + e) , and 
bounding each term 5 {s,s + e) by 5 (e), the diagonal portion of E ([X, m]^ {T))^ can be bounded 
above, in absolute value, by 

1 (™-l)/2 T 

— V c, / dt dsd^"" (e) = est ■ 8-^5^"" (e) . 

^ j^Q J2e Jt-2e 

The hypothesis on the univariate 6"^ implies that this converges to as e tends to 0. The case of 
t < 2e works as easily. 



Step 2: off diagonal. The ofi'-diagonal contribution is the sum for j = 0, • • • , (m — 1) /2 of the 
terms 

-2, 



j dt j ds5^^s,s + e)6^^t,t + e)e' is,tr-'^^ (22) 
As we will prove below, the dominant term turns out to be J[m-i)/2,OD\ deal with it now. 
Step 2.1: term J(m-i)/2,OD- Denoting c= |c(m-i)/2|i we have 



Jj,OD = e "Cj 



J, 



{m-l)/2,OD 



< 



T 

dt 

2e JO 



t-2e 



ds\Q' {s,t)\ . 



We estimate the integral, using the formula ()2ip and Fubini's theorem: 



T 

dt 

2e Jo 



t-2e 



ds {s,t)\ = I dt I ds 

J2e Jo 



t-2e 



s+e rt+e 



fi (dudv) 



T ct-2e /-s+e r-t+e 



< / dt 

he Jo 
(■T+£ fv-e 



ds 



s Jt 



{dudv) 



v=2e Ju=0 
T+e rv-e 



< 



\fi\ (dudv) 



fi\ (dudv) 

/min{u,t—2e) 
_ _ ^ _ =max(0,n— e) 

ds dt 



ds dt 



v=2e Ju=0 

T+e i-v-e 



t=v—e J s=u—e 



v=2e Ju=0 



lul (dudv) 
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Hence we have 

i-T+e i-v—e 

J{m-i)/2,OD<c6'^"'''^ie) / \fi\idudv) 

Jv=2e Ju=0 

< c6'^"'-^ (e) (OD) , 
which again converges to by hypothesis as e goes to 0. 

Step 2.2: other Jj^oD terms. Let now j < {m — 1) /2. Using Schwartz's inequahty for all but one 
of the m — 2j factors G in the expression (|22p for Jj,oD, which is allowed because m — 2j > 1 here, 
exploiting the bounds on the variance terms via the univariate function 6, we have 

|J,,od|<^-^^ f^dt ds\@^s,t)r-'^-'\Q^s,t)\ 

£ J2e JO 

f-T ft-2e 

< S^'^-^e) CjE-'^ dt ds\e'{s,t)\, 

J2£ Jo 

which is the same term we estimated in Step 2.1. This finishes the proof of the theorem. ■ 

A typical situation covered by the above theorem is that of the Riemann-Liouville fractional 
Brownian motion. This is the process B^'^^ defined by B^'^^ (t) = (t — s)^"^^"^ dW (s). Its 
canonical metric is not homogeneous, but we do have, when H € (0, 1/2), 

\t-s\^ <5{s,t) <2\t-s\^ , 

which implies, incidentally, that B^'^^ has the same regularity properties as fractional Brownian 
motion. To apply the theorem, we must choose H > l/(2m) for the condition on the variances. 
For the other condition, we calculate that fi (dsdt) = 2H (1 — 2H) \t — s]"^^^^ dsdt, and therefore 

/i (OD) = (OD) = CH s^^-'^ds < CHTe^^-\ 

This quantity is bounded above by e^i+^/(^"^) as soon as H > 1/ (2m), of course, so the strict 
inequality is sufficient to apply the theorem and conclude that B^'^^ then has zero mth variation. 

One can generalize this example to any Gaussian process with a Volterra-convolution kernel: 
let 6'^ be a univariate increasing concave function, differentiable everywhere except possibly at 0, 
and define 

X (t) = [^j (r) dW (r) . 

Then we can easily show that the canonical metric 5"^ is,t) of X is bounded above by 26"^ {\t — s\), 
and also below by 6'^ {\t — s\), and similar calculations to the above show that X has zero mth 
variation as soon as (5^ (r) = o (r^/^-^'")) . Hence there are inhomogeneous processes that are more 
irregular than fractional Brownian for any H > 1/ (2m) which still have zero mth variation: use 
for instance the X above with 6^ (r) = r^/C^"^) / log (1/r). 

3.4 Non-odd powers 

When m > 1 is not an odd integer, recall that to define the mth odd variation, we use the convention 
(x)"^ = Ixl"^ sgn (x), which is an odd function. The idea here is to use the Taylor expansion for this 
function up to order [m], with a remainder of order [m] + 1; it can be expressed as the following 
lemma. 
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Lemma 8 Fix m > 1 and two reals a and b such that \a\ > Let (™) denote the formal binomial 
coefficient m [m — 1) • • • (m — k + 1) / {k [k — 1) ■ ■ ■ 1) and let {x) := |a;| sgn{x). Then, for all reals 
a,b, 

1. if\b/a\ < I, 

H-i 



I — n ^ ^ 



\ I \2m—k ik I I 1 2m ^ 

a) \a\ b + \a\ frn-i 



k=0 

2. and if \a/b\ < 1 



k=o ^ ^ 



a 

m 



67' 



where for all \x\ < 1, \ fm {x)\ < Cm \x\^"^'^^^ where Cm depends only on m. 

This lemma is proved in the Appendix. When we apply it to the question of finding Gaussian 
fields with zero odd mth variation, we are able to prove that the sufficient condition of Theorem [6] 
still works. The result is the following. 

Theorem 9 Let X be as in Theorem\^ and m be any real number > 1. Then d {r) = o (r^/^^™)) 
implies X has zero odd mth variation. 

Proof. Step 0: setup. Recall that with Y = X {t + e) - X {t) and Z = X (s + e) - X (s), we 
have 



E 



{[XMe{T)f]=^^ f fdtdsE[{Yr{Zr 
J ^ Jo Jo 



Now introduce the shorthand notation cj^ = Var [Y], = Var [Z], and 9 = E [XY] = 0. Thus 
Y = aM where M is a standard normal r.v.. We can write the "linear regression" expansion of Z 
w.r.t. Y, using another standard normal r.v. N independent of M: 

9 

Z = -M + rrN 
a 

where 

Note that r is always well-defined and positive by Schwartz's inequality. Therefore 

/ 9 \ 

{Ynzr = a"" {Mr i^M + rrNj = sgn (0) a^'" j^p'" (a)"^ (a + 6)' 



where 



a := —M and b := rrN. 
a 
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Applying Lemma [8l we get that (Y)"^ (Z)™" is the sum of the following four expressions: 

[m]-l 



A 



fc=0 



sgn 



k+l 



(0) 



\m-k ^k^k^k 



A' := sgn {6)l\rra 



B := 1, 



fc=0 



ni 



|,.|2m. frTaN\ 



sgn 



fc+2 



(^) 



m+fc i^im— fc 



B' :=SKn(g)l|,...ivK Jmrr {MN)'^ U 



OM 



(23) 
(24) 
(25) 
(26) 



S'iep i; cancellations in expectation calculation for A and B. In evaluating the e-mth variation 
E ([X, m]e (T))^ , terms in A and B containing odd powers of M and N will cancel, because of 
the symmetry of the normal law, of the fact that the indicator functions in the expressions for A 
and B above are even functions of M and of N, and of their independence. Hence we can perform 



the following calculations, where a^n^k '■= E 
positive constants depending only on m and k. 



|^|2m-fc|^|fc 



and h 



'm,k 



E 



\Mr+''\N\ 



m—k 



are 



Step 1.1: expectation of A. In this case, because of the term N^, the expectation of all the terms 
in (j23l) with k odd drop out. We can expand the term using the binomial formula, and then 
perform a change of variables. We then have, with n = [m] — 2 when [m] is even, or n = [m] — 1 
when \rn\ is odd. 



IE MIS E T l« 

u—n v / 



m~k kk^ 

a T a, 



k=0 

k even 



k/2 



n/2 n 

^|^r-2^'(ar)2^^ V 



k=2j 
k even 



k 
k 



k/2-j 



(_l)'=/2-i 
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(^m,k 



n/2 

sEi 

i=o 



im—2j 



(ar) 



2j 



where Cm,j are positive constant depending only on m and j. In all cases, the portion of E ([X,m]e {T)y 
corresponding to A can be treated using the same method as in the proof of Theorem [6l . This 
lowest power is attained when j = n/2, i.e., when [m] is even we have |0|2+™'~M^ g^j^^j when [m] is 
odd, we have 1^1^+™-"!™!. In both cases, the power is greater than 1. All other values of j correspond 
of course to higher powers of \6\. This means we can use Schwartz's inequality to get the bound, 
valid for all j, 

and we are back to the situation solved in the proof of Theorem[6j Thus the portion of E {[X, mj^ (T))^ 



o r 



l/{2m)\ 



corresponding to A tends to as e ^ 0, as long as 5 (r) 

Step 1.2: expectation of B. This portion is dealt with similarly. Because of the term sgn'^"'"^ 
the expectation of all the terms in (j25p with k even drop out. Contrary to the case of A, we do not 
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need to expand r*" in a binomial series. Since k is now > 1, we simply use Schwartz's inequality 



to write \6\'' < \9\ {<tt)^ ^ . Of course, we also have r < 1. Hence 



|E[S]| 



k=l ^ ^ 

k odd 

<\e\{arr~' E ^-'HT • 
k=i ^ ^ 



<J T 



(27) 



We see here in all cases that we are exactly in the same situation of the proof of Theorem [6] (power 
of 1^1 is 1^1^). Thus the portion of E ([X, m]^ (T))'^ corresponding to B converges to as soon as 

Step 2. The error term A' . For A' given in ()24p . we immediately have 

/ rraN 



|E [A'] 



E 



l|rz^|<ll^n^l'"/m-l 



1 1 rra 



<cm\orB 

= c^l0r"t™Urrf7)[™]E 



2m 



rraN 


[m] 


9M 





1| rra-N I ^-t \M 



2m- H |jY|H 



We see here that we cannot ignore the indicator function inside the expectation, because if we did 
we would be left with |^| to the power m — [m], which is less than 1, and therefore does not allow 
us to use the proof of Theorem O 

To estimate the expectation, let x = ^j^. We can use Holder's inequality for a conjugate pair 
p, q with p very large, to write 



E 



1,^1^1 IMI^"^-!™] lA^jt™]] < P^/^lxNl < \M\] E 



,1/p 



|M 



2mp-[m]p [m]p 



The second factor in the right-hand side above is a constant Cm,p depending only on m and p. For 
the first factor, we can use the following standard estimate for all y > 0: e~^'^^^dz < y~^e~y^ ^'^ . 
Therefore, 



V[\xN\ < \M\] = 2 



~ [|M| > ™| < , /I e-«v. + ,/I r <,„£rf^!, 

V27r V vr Jo V vr Jy^ ux 

-1,2/2 



-I roo -| / 1 1 poo 1 

< ^27^- + / dn— e-"'"'/2 ^ - + - / di^-e 



c 

X 



where c is a universal constant. 

Now choose q so that m — [m] + 1/q = 1, i.e. q = (1 — m + [m]) which exceeds 1 as long as 
m is not an integer. Then we get 



|E [A'] I < Cm |er-["^l (rTa)H cm,p 



m,p 1 = CmCm,pC^/'*' |6'| (rrcj) 
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and we are again back to the usual computations. The case of m integer is dealt with in Step 4 
below. 



Step 3. The error term B' . For B' in ()26p . we have 



|E [B'] \ < c^c7™T'"r"E 



iMiVl 



9M 



rroN 



\m]+l 



The expectation above is a constant depending only on m as long as m is not an integer. The 
case of m integer is trivial since then we have B' = 0. Now we can use Schwartz's inequality to say 
that < (err) ['"I, yielding 

\m~l I /II 1 



\E[B']\ < c„ 



c' r''' 



{ary 



This is again identical to the terms we have already dealt with, but for the presence of the negative 
power on r. We will handle this complication by showing that r can be bounded below by a 
universal constant. First note that integration on the e-diagonal can be handled by using the same 
argument as in Steps 1 and 2 of the proof of Theorem [6l Thus we can assume that t > s + 2e. Now 
that we are off the diagonal, note that we can write without loss of generality, using the mean value 
theorem, that 6 = (5^) (^) for some ^ in [t — s — e,t — s + e\: here ((5^) is the second derivative 
of the convex function 5^. Now because we are assuming that our process X is more irregular 

than Brownian motion, we can assume without loss of generality that 



Indeed, the case of 5"^ ir) 



.1/2 



(52)" (r) <(3/4)r-252 



6' ir) 



„2H 



(fBm with H = 1/4) yields an equality, any other function 



with H < 1/2 yields a smaller proportionality constant than 3/4, and beyond the 



Holder scale, we actually have 



o r 



-2x2 



Also assuming 



[5^ 



IS mcreasmg 



without loss of generality, since ar = 5'^ (e), we can now say that off the diagonal, \6/ ((Tr)| < 
(<^^)" (^) (^) ^ 3/4, which proves that r bounded below uniformly (off the e-diagonal) by 
(1 — 9/16)^^^. Hence, in the inequality for |E [B']\ above, the term y-^-M-i can be absorbed into 
the m-dependent constants, proving the corresponding contribution to E {[X, m]^ (T))^ converges 
to as soon as (5^ (r) = o (r^/^^™-)). 

Step 4- The case of m integer. Of course, we already proved the theorem in the case m odd. Now 
assume m is an even integer. In this special case, we do not need to use a Taylor expansion, since the 
binomial formula has no remainder. Moreover, on the event |^-ffp| < 1, sgn (a + 6) = sgn (a). Thus 
A' = with the understanding that we replace A by the full sum for k = thru m. Recalculating 
this A we get 



A = cT^sgn {9) sgn (M) iMp ^ Ij^'^N 



k=0 



rm—k 



-I k 



a 



rr 



-im—k 



k=a ^ ^ 



{M)\M\^^''N 



m+k ]\Tm—knk„m—k^m—k^m—k 



a 



T 



Here, when we take the expectation E, all terms vanish since we have odd functions of M for k 
even thanks to the term sgn^"*""^ (M), and odd functions of N for k even thanks to the term ]^^~^ , 
I.e. the term corresponding to A is entirely null when m is even. ■ 
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4 Non- Gaussian case 



Now assume that X is given by ([5]) and M is a square-integrable (non-Gaussian) martingale, m is 
an odd integer, and define a positive non-random measure /i for s = (si, S2, • • • , Sm) £ [0, T]^ by 

fi (ds) = fi {dsids2 ■ ■ ■ dsm) = E [d [M] (si) d [M] (ss) ■■■d [M] , (28) 

wliere [M] is tlie quadratic variation process of M. We make the fohowing assumption on fj,. 

(M) The non-negative measure /i is absolutely continuous with respect to the Lebesgue measure ds 
on [0, r]™ and K (s) := dfi/ds is bounded by a tensor-power function: < (si, S2, • • • , Sm) < 
(si) (§2) • (sm) for some non-negative function T on [0, T]. 

A large class of processes satisfying (M) is the case where M (t) = /q H (s) dW (s) where H E 
-L^ ([0, T] X Q) and PV^ is a standard Wiener process, and we assume E [ff^*" (t)] is finite for all t E 
[0, T]. Indeed then by Holder's inequality, since we can take (s) = E [H'^ (si) (^2) • • • (sm)] , 

we see that T (s) = (E [H^"" (t)])^/^^""^ works. 

We will show that the sufficient conditions for zero odd variation in the Gaussian cases generalize 
to the case of condition (M), by associating X with a Gaussian process. We let 

G{t,s) = T{s)G{t,s) 

and ^ 

Z{t)= [ G {t, s) dW (s) . (29) 
Jo 

We have the following. 

Theorem 10 Let m be an odd integer > 3. Assume Z defined in 1^29^) satisfies Condition (M) and 
the hypotheses of Theorem or Theorem relative to a univariate function 5. Assume that for 
some constant c > 0, and every small e > 0, 

/ dt I ds I AGt (u) AGs (u) du < ce5^ (2e) . (30) 

Jt=2£ Js=0 Ju=0 

Then X has zero m.th variation. 
Proof. 

Step 0: setup. We use an expansion for powers of martingales written explicitly on page 14 of 
[lOj . For Y a fixed square-integrable martingale, for any integer k E [0, [m/2]], with the set of 
permutations oi m — k defined as those for which the first k terms are chosen arbitrarily and the 
next m — 2k terms are taken as the complement of the first k terms chosen, we define the process 
Y^^£ (denoted in the above reference by a^) by setting, for each a E and each ^ = 1, 2, • • • ,m — k, 



We then have for all t E [0, T] 



[Y]{t) if a(^)E{l,2,--- ,k} 
Y{t) [i a{l) £ {k + 1,--- ,m- k} . 



{Ytr = E ^ E / / • • • / ^^--1 ("1) ^^-.2 (^2) • • • dY^^^^k {um-k) 



k=0 o-eE 
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We use this formula to evaluate 



[X, m]e (T) = - r ds{X{s + e)-X {s))' 
^ Jo 



by noting that the increment X {s + e) — X (s) is the value at time T of the martingale Yt 
Jq AGs {u) dM (n) where we set 



AGs {u) ■.= G{s + £,u) -G{s., 



u) 



Hence 

{X{s + e)-X{s)r 



EwE/ / •••/ ^^[M]K(i))|AG,(«,(i))|^..d[M](n,(,))|AG,K(fc))|- 

dM {u„(k+l)) AGs (^^<7(fc+l)) ■■■dM {u„(rn-k)) {Ua(rn-k)) ■ 



Therefore we can write 



= - E ^ E / / ■ ■ ■ / ■■■dW\ [u^ik+i)) ■■■dM [u^^m-k)) 

^ k=0 ^ .gE^ 

[AG. (n^(fc+i)) ; • • • ; AG. (M^(^_fc)) ; AG. (u^ci)) ; {^^(1)) r ■ ■ ; AG. (^^.(fc)) ; AG. (ua(fe))] • 



where we used the notation 



[fl, /2, • • • , fm] ■■= / /l (s) h («)••• /m (s) ds. 

Jo 



To calculate the expected square of the above, we will bound it above by the sum over k and a 
of the expected square of each term. Writing squares of Ricmann integrals as double integrals, 
and using Ito's formula, each term's expected square is thus, up to (m, A;)-dependent multiplicative 
constants, equal to the expression 

/ / / ■■■/ / 

^ Jum-k=0 Ju'^_i,=0 Ju,n-k-i=0 Ju'^_^_-^=0 Jui=oJu[=0 
B[d[Mf' ,n^ik))d[Mf' (<(!)'••• ^<ik))d[Mf^"''''^ ,n^im-k)) 

■ [AG. (n^(fc+i)) ; • • • ; AG. {u^(m-k)) ; AG. (u^d)) ; AG. {u„(i)) ;■■■ ; AG. {u^(k)) ; AG. {u^(k))] 

■ AG. {u,^k+i)) ; • • • ; AG. {u,^m-k)) ; AG. ; AG. ; • • • ; AG. ; AG. 

(31) 

modulo the fact that one may remove the integrals with respect to those u^s that are not rep- 
resented among {u'^^-^y--- ,m^j^^,j}. The theorem will now be proved if we can show that for all 

k e {0, 1, 2, • • • , [m/2]} and all a € S^, the above expression K = i^^n^fe^o-tends to as e tends to 
0. 
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A final note about notation. The bracket notation in the last two lines of the expression (|3ip 
above means that we have the product of two separate Riemann integrals over s £ [0,T]. Below 
we will denote these integrals as being with respect to s G [0,T] and t £ [0, T]. 

Step 1: diagonal. As in Steps 1 of the proofs of Theorems [6] and [3 we can use brutal applications 
of Schwartz's inequality to deal with the portion of Km,k,a in (|3ip where |s — t| < 2e. The details 
are omitted. 

Step 2: term for k = 0. When A; = 0, there is only one permutation a = Id, and we have, using 
hypothesis (M) 



K, 



m,0,ld 9 



< 



- / ••• / E[d[M]^™(ni,... ,n„)] •[AG.(ni);--- ; AG. (n^)r 

£ Jum=0 Jum-1=0 Ju\=0 
T 



«m-fc=0 Ju. 
T 



m-1—0 
Um — k 



^m — k 



=0 



U2 

ui=0 

Ml=0 



(ni) {u2) • • • (um) [AG. (m) ; • • • ; AG. (um)] duidu2 ■ ■ ■ dur. 



AG. (ni);-- - ; AG. (n„) 



duidu2 ■ ■ ■ dUr^ 



This is precisely the expression one gets for the term corresponding to A; = when M = W, i.e. 
when X is the Gaussian process Z with kernel G. Hence our hypotheses from the previous two 
theorems guarantee that this expression tends to 0. 

Step 3: term for k = 1. Again, in this case, there is only one possible permutation, a = Id, and 
we thus have, using hypothesis (M), 
rT 



KmAJd = 



Mm-1=0 Jum~2=0 



E 



ui=0 Jm'i=0 



d [M] (m) d [M] (n'l) d [M]®("^-2) {U2, • • • , n„_i) 



• [AG. (U2) ; • • • ; AG. {um-i) ; AG. (m) ; AG. (m)] • [AG. (na) ; • • • ; AG. ; AG. (n'l) ; AG. {u'^)] 



< 



U2 /•U2 



Ul=0 Ju[=0 



duidu'idu2 ■ ■ ■ dum-iT'^ (ui) (n'l) {u2) • • • (um) 



• [|AG|. {U2) ; • • • ; |AG|. (um-i) ; |AG|. {m) ; |AG|. {m)] • [|AG|. {u2) ; • • • ; |AG|. (n„_i) ; |AG|. {u[) ; |AG|. «)] 



1 

AG 



T 



Um~l=0 Ju,n~-2=0 



duidu'idu2 ■ ■ ■ du. 



1 



ui=0 Ju'=Q 



AG 


{u2) ; • • • ; 


AG 


(um-i) ; 


AG 


^ (^^i) ; 


AG 





im) ; • • • ; 



AG 



(um-i) ; 



AG 



AG 



i) 



Note now that the product of two bracket operators [•••][•••] means we integrate over < s < t—2e 
and 2e <t <T, and get an additional factor of 2, since the diagonal term was dealt with in Step 1. 

In order to exploit the additional hypothesis (j30p in our theorem, our first move is to use Fubini 
by bringing the integrals over ui all the way inside. We get 

f-T 



T 



«m-l=0 J Um-2=0 

t-2e 

ds dt 



U3 



duo ■ ■ ■ du 



m—l 



U2=0 



t=2e Js=Q 
U2 fU2 



AGs{u2) ■■■ AG, AGt(n2) ••• AGt{um-i) 

(■^2 (■U2 , _ X 2 / _ \ 2 

/ / duidu'JAG,{u{)) (AGt(n; 
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The term in the last hne above is trivially bounded above by 

cT cT 



[ [ duidu'^ (aG, (m))' (/\Gt (n'l 



precisely equal to Var [Z {s + e) — Z (s)] Var [Z {t + e) — Z (t)], which by hypothesis is bounded 
above by 5^ (e). Consequently, we get 



iriAJd 



< 2 



Um-l 



U3 



P II I 

^ JUm-l=0 JUm,-2=0 JU2=0 

T rt-2e 



dU2 ■ ■ ■ dUm-1 



ds dt 



AG, (na) ••• AG, (u^-i) AG* (ua) ••• AGt{um-i 



lt=2£ Js=0 

We get an upper bound by integrating all the tij's over their entire range [0,T]. I.e. we have, 
6^ (e) 



m,l,Id 



< 



dt / ds 

t=2e Js=0 



T rT 



JO 
T 



du^ ■ ■ ■ duj 



m—1 







AG,(n3) ••• AG,(u„_i) AGtiu-s) 



AGt (Um^i] 



U2=0 



AGt{u2) AGs{u2 

t-2e 



dt 

Jt=2£ Js=o 



ds 



dU2 

T 



du 







AG, (u) AGt{u 



m— 3 



«3 



«2=0 



AGt{u2) AGs{u2 



duo 



Now we use a simple Schwartz inequality for the integral over u, but not for U2- Recognizing that 
AGs (u) du is the variance Var [Z (s + e) — Z (s)] < 5^ (e), we have 



Km,l,Id < 2 

< 2 



6' (e) 






r 

/t=2e 


^4+2m- 


-He) 



dt 



t-2e 



ds 



T 



du 







AG, [u] 



2\ '"-^ 



lt2=0 



T 



T 



/■t~2e 
dt / ds 
t=2e Js=0 Ju2=0 



AGt{u2) AGs{u2) 



AGt{u2) AGs{u2) 

dU2- 



duo 



Condition (j30p implies immediately K^^ijd ^ i^^™ (2^) £ ^ which tends to with e by hypothesis. 

S'tep 4- k > 2. This step proceeds using the same technique as Step 3. Fix k > 2. Now for each 
given permutation a, there are k pairs of parameters of the type {u,u'). Each of these contributes 
precisely a term (e), as in the previous step, i.e. 6^^ (e) altogether. In other words, for every 
a E S^, and deleting the diagonal term, we have 



< 2^4^ I dt I ds 



t=2e Js=0 



Jo 



'■k+2 



dUk+l ■ ■ ■ dUm-k 



ds 



AGs (uk+i) 



AGs (Ur^ 



Since k < (m — 1) /2, there is at least one integral, the one in u^+i, above. We treat all the re- 
maining integrals, if any, over Uk+2, • • • jUm-k with Schwartz's inequality as in Step 3, yielding a 



26 



contribution 52(m-2fc-i) ^^^^ rpj^^ 

remaining integral over Uk+i yields, by Condition (|30p . a contri- 
bution of (5^ (2e) e. Hence the contribution of Km,k,a is again (2e) e~^, which tends to with e 
by hypothesis, concluding the proof of the Theorem. ■ 

We state and prove the next proposition, in order to illustrate the range of applicability of 
Theorem [TOl The presence of the indicator function ls<f in the expression for G below is typical 
of most models, since it coincides with asking that Z be adapted to the filtrations of W , which is 
equivalent to X being adapted to the filtration of M. In the case of irregular processes, which is 
the focus of this paper, the presence of the indicator function makes G non-monotone in both s 
and t, which creates technical difficulties. Examples of non-adapted irregular processes are easier 
to treat, since it is possible to require that G be monotone. We do not consider such non-adapted 
processes further. More specific examples of adapted processes are given after the proof of the next 
proposition. 

Proposition 11 Let X he defined by 1^ and assume condition (M) holds. Let X (t) = jT T (s) G {t, s) dW {s), 



where W is a Wiener process. Assume that G [s] 
for all s, t, 

G {t, s) = ls<t g (t, s) = ls<t \t 



T (s) G {t, s) can be bounded above as follows: 



in which the bivariate function f (t, s) is positive and bounded as 

\f{t,s)\<f{\t-s\) 

where the univariate function f (r) is increasing, and concave on 
where g has a second mixed derivative such that 



with lim^^o / ('^) = 0, and 



dg_ 
dt 



it,s) 



+ 



dg_ 

ds 

d'g 



dsdt 



{t,s) 



{t,s] 



<c t 



< c t 



|l/(2m)~3/2 



|l/(2m)-5/2 



Also assume that g is decreasing in t and the bivariate f is increasing in t. Then X has zero cubic 
variation. 

Proof. Below the value 1/ (2m) — 1/2 is denoted by a. We now show that we can apply 
Theorem [7] directly to which, by Theorem IIOI is sufficient, together with Condition (|30p . to 
obtain our desired conclusion. We will prove Condition (|18p holds in Step 1; Step 2 will show 
Condition ([19]) holds; Condition ([30|) will be established in Step 3. 

Step 1. Variance calculation. We need only to show 5^(s,s + e) = o (e^/™) uniformly in s. We 
have, for given s and t = s + e 



6"^ (s,s + e) 



dr 



G{t,r)-G{s,r 

r \{s + e-rr f{s + e,r 
Jo 



) — (s — r)" / {s, r)\ dr 



+ 



S+£ 



\s + e 



1 2a 1-2 



f {s + e,r)dr 



(32) 



--:A + B. 
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Since (s + e, r) < / (s + e — r) and the univariate / increases, in B we can bound this last 
quantity by / (e), and we get 



B < f (e) [ r^^dr = Sf (e) e 
Jo 



The term A is sHghtly more dehcate to estimate. By the fact that / is increasing and g is 
decreasing in t, 

A< r f{s + e,r)\{s + e-r)'' -{s-r)''fdr= f f {e + r)\r^ - {r + eT\^ dr 
Jo Jo 

f {e + r) |r" - (r + e)"|^ dr + [ f (e + r) |r" - (r + ef\'^ dr 



=: A1+A2. 

We have, again from the univariate /'s increasingness, and the hmit hm,._+o / (^) = 0, 

Ai < f (2e) |r° - (r + e)"!^ dr = est f {2s) e^^+i = o (e^/™) . 

For the other part of A, we need to use /'s concavity at the point 2e in the interval [0, e + r] (since 
e + r > 2e in this case), which implies f (e + r) < f (2e) (e + r) / (2e). Also using the mean- value 
theorem for the difference of negative cubes, we get 



/s ps 
f {e + r) r'^°'~^dr < est ef (2e) J (e + r) r 

< est ef {2e) J' r^^-^ = est e^^+i/ {2e) = o (e^' 



This finishes the proof of Condition (jlSp . 

Step 2. Covarianee calculation. We first calculate the second mixed derivative d'^o^ /dsdt, where 5 
is the canonical metric of Z, because we must show (OD) < e^", which is condition (jl9p . and 
fi (dsdt) = ds dt dH"^ /dsdt. We have, for < s < t — e, 

i9{t,s - r) - g{s,s - r)f dr + g'^{t,r)dr 
Jo Js 

=:A + B. 

We calculate 

|^(t,s) = 2^(t,0) {g{t,0)-g{s,0)) 

f'\dg . . / dg . . dg , . dg . 

2, 



r 5 5 

+ y 2 (5 (t, s - r) - 5 (s, s - r)) (t, s - r) dr. 



'0 

A1 + A2 + A3 
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and 



d^B , , , ^ do , ^ 



Next, we immediately get, for the portion of iOD) corresponding to B, using the hypotheses 
of our proposition, 



T 



dt 



t-e 



ds 



dsdt 



it,s] 



t-e 

<2c I dt I dsf{\t-s\)\t-s\"'\t-s 

■T 



a-1 



< 2c 



dt = est e^", 



which is of the correct order for Condition ()19p . For the term corresponding to Ai, using our 
hypotheses, we have 



dt 



ds \Ai\ < 2 



dt 



t-e 



ds f 



where S^t,s is in the interval (s,t). Our hypothesis thus implies ^ {^t,s,^) 



t - s\ 
< s", and hence 



dt 



t-e 



T 



t-e 



ds \Ai\ <2T I dt I ds s"-^t"-^ = 2Ta-^ 



T 



dt {t - e)" < a-'T 



-2rpl+2a 



This is much smaller than the right-hand side e^" of Condition ()19p . since 2a = 1/m — 1 < 0. The 
terms A2 and ^3 are treated similarly, thanks to our hypotheses. 

Step 3: proving Condition (fgOj). In fact, we modify the proof of Theorem [TOl in particular Steps 3 
and 4, so that we only need to prove 



dt 



t=2e Js=0 



t-2e 



ds 



u=0 



AGtiu) AGs{u] 



(33) 



instead of Condition (j30p . Indeed, for instance in Step 3, this new condition yields a final contri- 
bution of order (52m-2 ^^-j ^-2gi/m+i^ With the assumption on 5 that we have, 5 {e) = (e-^/^^™)), 
and hence the final contribution is of order o (e(2r?i-2)/(2rra)-i+i/m'j _ This proves that the 

conclusion of Theorem [TU] holds if we assume (|33p instead of Condition (|3Up . 
We now prove (j33p . We can write 



dt 



t-2e 



ds 



t=2e J s=Q 



dt 



u=0 



AG. iu] 



AG, {u] 



du 



t=2e Js=Q 



t-2e 



ds 



+ 



dt 

t=2e Js=0 

A + B. 



t-2e 



ds 



\g{t + £,u) - g{t,u)\ \g{s + e,u) - g{s,u)\du 

he 

\g{t + e,u) - g{t,u)\ \g{s + e,u)\du 



For A, we use the hypotheses of this proposition: for the last factor in A, we exploit the fact 
that g is decreasing in t while / is increasing in t; for the other factor in A, use the bound on dg/dt; 
thus we have 



A < 



dt 



t=2e Js=0 



t-2e 



£\t- s|"-' ds f{s + e, u) {{s -u)" -{s + e- n)") du. 
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We separate the integral in u into two pieces, for n G [0,s — e] and u G [s — e,s]. For the first 
integral in n, since / is bounded, we have 



s—e fs—e 

""1 ^o, <^ II #l| ^ A+a 



/ f{s + e,u) {{s -ur-{s + e- uf) du<\\f\\^e {s - uf^' du < \\f\\^ c.e 
'0 Jo 



For the second integral in ti, we use the fact that s — u + e > e and s — u < e implies s — u + e > 
2{s — u), so that the negative part of the integral can be ignored, and thus 

/ (s + e, u) {{s -ur-{s + e- n)°) du < T {s - uf du = c„el+^ 



s—e 



which is the same upper bound as for the other part of the integral in u. Thus 

A < est £2+" / dt \t- s\"-^ ds < est £2+" / dte" < est e^+a" = est e^/™+\ 

Jt=2e Js=0 Jt=2e 

which is the conclusion we needed at least for A. 

Lastly, we estimate B. We use the fact that / is bounded, and thus \g {s + £,u)\ < \\f\\^ \s + e — u\", 
as well as the estimate on the derivative of g as we did in the calculation of A, yielding 

B<\\f\\^e dt ds|t-s-er~M \s + e - u\'^ du 

Jt=2e Js=Q J s 

este"^^ I dt / ds\t-s-er-^ 



lt=2e Js=0 

pT i-t-2e 

< 2i+l"lcst e"+2 / dt ds\t- s|"-^ < est e2°+2 = est e'^/"'+\ 

Jt=2e Js=0 

This is the conclusion we needed for 5, which finishes the proof of the proposition. ■ 

The above proposition covers a very wide variety of martingale-based models, in the sense that 
they can be quite far from Gaussian models. We describe one easily constructed class. Assume 



that M is a martingale such that E 



\d[M] /dt\^'^ 



is bounded above by a constant c uniformly 



in i < T. Let G{t,s) = GnLfBm {t,s) := ls<t \t — s|^/(2™')~^/2+" for some a > 0; in other words, 
G is the Brownian representation kernel of the Riemann-Liouville fractional Brownian motion with 
parameter H = l/(2m) — a > 1/ (2m). Or more generally, assume that G is bounded above by 
GRLfBm-, and assume the two partial derivatives of G, and the mixed second order derivative of G, 
are bounded by the corresponding derivatives of GjiLfBm- Then by the previous proposition, X 
defined by ([5]) has zero mth variation. 

The most quantitatively significant condition, that the univariate function 5 (e) corresponding 
to G be equal to o (e^^^^™^), can be interpreted as a regularity condition. In the Gaussian case, 

it means that there is a function /(e) = o ^e^/^^^Mog^^^ (1/e)^ such that / is an almost-sure 

uniform modulus of continuity for X. In non-Gaussian cases, similar interpretations can be given 
for the regularity of X itself, provided enough moments of X exist. If X has fractional exponential 
moments, then the function / above will also serve as an almost-sure uniform modulus of continuity 
for X, provided the logarithmic correction term in / is raised to an appropriately higher power 
than 1/2. This can be established by using the non-Gaussian regularity theory in [23j. If X has 
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standard moments of all orders, then one can replace /(e) by for any q > 0. This 

is easily achieved using Kolmogorov's continuity criterion. If X only has finitely many moments, 
Kolmogorov's continuity criterion can only guarantee that one may take a greater than some > 0- 
We do not delve into the details of these regularity issues in the non-Gaussian martingale case. 



5 Stochastic calculus 

In this section, we investigate the possibility of defining the so-called symmetric stochastic integral 
and its associated Ito formula for processes which are not fractional Brownian motion; fBm was 
treated in |12]. We concentrate on Gaussian processes under hypotheses similar to those used in 
Section 13.31 (Theorem [7]). 

The basic strategy is to use the results of [l2]. Let X be a stochastic process on [0, 1]. According 
to Sections 3 and 4 in [12] (specifically, according to the proof of part 1 of Theorem 4.5 therein), if 
for every bounded measurable function g on R, and every integer n > 1, the limits 



lim i du (X„^., - X^r^' g ( ^"+-+^" ) = 



(34) 



hold in probability, then for every t G [0,1] and every / G (R), the symmetric ("generalized 
Stratonovich" ) stochastic integral 



^* /' (X„) X (odu) =: lim ^ du (X,+, - XJ /' (^ ^-+^+^- ^ = q 



(35) 



exists and we have the Ito formula 



/ (X) = / (Xo) + /* /' (X„) X (odu) . (36) 
Jo 

Our goal is thus to prove ()34p for a wide class of Gaussian processes X, which will in turn imply 
the existence of ([35]) and the Ito formula ([36]) . 

If X has homogeneous increments in the sense of Section [3.2| meaning that E (X^ — Xj)^ = 

6^ (i — s) for some univariate canonical metric function 6, then by using g = 1 and our Theorem [6l 
we see that for p4p to hold, we must have 5 (r) = o (r^/^) . In order to go beyond the homogeneous 
case, we notice that (I34p for g = 1 is the result of our non-homogeneous Theorem[71 so it is necessary 
to use that theorem's hypotheses, which include the non-homogeneous version of 5(r) = o(rV6). 
But we will also need some non-degeneracy conditions in order to apply the quartic linear regression 
method of [12]. These are Conditions (i) and (ii) in the next Theorem. Condition (iii) therein is 
essentially a consequence of the condition that 6'^ be increasing and concave. These conditions are 
all further discussed below. 



Theorem 12 Let m > 3 be an odd integer. Let X be a Gaussian process on [0, 1] satisfying the 
hypotheses of Theorem This means in particular that we denote as usual its canonical metric 
by 5^ {s,t), and that there exists a univariate increasing and concave function 6'^ such that 5 (r) = 
o (r-^/*^^™)) and 5"^ is,t) < 5"^ {\t — s\). Assume that for u < v, the functions u Var [X^] =: Qu, 
V ^ 5"^ {u,v), and u ^ —5"^ iu,v) are increasing and concave. Assume there exist constants c,c' 
such that 
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(i) < Q 



(ii) (2 + c)QuQ (n, v) + {l- c) QuQv + (2 - c) > q2 ^) ^ 

(iii) i/iere exist a > 1 and 6 < 1/2 sitc/i that for all u G [e, 1], 

5 iau) — 5 (u) , (5 (u) 



(a — 1) n 



u 



(37) 



r/ien /or every bounded measurable function g on R, 



lim 

£^0 



As mentioned above, this theorem immediately imphes the fohowing. 



Corollary 13 Assume the hypotheses of Theorem with m = 3. We have existence of the 
symmetric integral in ( f35|] . and its ltd formula h3(j\) . for every f ^ (R) and t G [0, 1] . 

Before proceeding to the proof of this theorem, we discuss its hypotheses. We refer to the 
description at the end of Section 13.31 for examples satisfying the hypotheses of Theorem [71 these 
examples also satisfy the monotonicity and convexity conditions in the above theorem. 

Condition (i) is quite weak, since c' can be chosen arbitrarily close to 0. The hypotheses of 
Theorem [7] imply that Qu < S"^ {u) , and Condition (i) simply adds that these two quantities should 
be commensurate, which is a weak non-degeneracy assumption, as opposed, say, to requiring that 
c' not be too small, which would be a stronger coercivity assumption. 

Condition (ii), while ominous looking, is not much more restrictive than Condition (i). Indeed, 
by brutally applying Schwartz's inequality leading to Q{u,v)^ < QuQv, Condition (ii) becomes a 
trivial consequence of 

cQv < (2 + c)Q {u, v) + {2- c) Qu. 

If the covariance function Q {u, v) is everywhere non-negative, this last inequality is then implied 
by the existence of a constant c such that cQy < Qu, which is equivalent to asking that the 
ratio Qu/Qv be bounded above and below by positive constants. Similarly to Condition (i), this 
commensurability is a weak non-degeneracy assumption. The non-negativity of Q {u, v) holds for 
many processes: since X (t) = £ G (t, s) dW (s) is defined via a Volterra representation w.r.t. 
Brownian motion, if the kernel G is non-negative, then so is Q {u,v). 

Lastly, Condition (iii) represents a strengthened concavity condition on the univariate function 
6. Indeed, the left-hand side in (f37|) is the slope of the secant of the graph of 6 between the points 
u and au, while the right-hand side is b times the slope of the secant from to u. If b were allowed 
to be 1, (iii) would simply be a consequence of convexity. Here taking 6 < 1/2 means that we are 
exploiting the concavity of (5^; the fact that condition (iii) requires slightly more, namely b strictly 
less than 1/2, allows us to work similarly to the scale S (r) = with H < 1/2, as opposed to 
simply asking H < 1/2. Since the point of the Theorem is to allow continuity moduli which are 
arbitrarily close to r^^^, Condition (iii) is hardly a restriction. 

Proof of Theorem I12L 
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Step 0: setup. The expectation to be evaluated is written, as usual, as a double integral over 
{u,v) G [0, 1]^. For e > fixed, we define the "off-diagonal" set 

De = {{u,v) £ [0, 1]2 : e^-P <u< v - e^-P < f < l} . 

Using the boundedness of g and Schwartz's inequality, thanks to the hypothesis 5 {r) = o (r^/^^™)), 
the term corresponding to the diagonal part (integral over W^) can be treated identically to what 
was done in p2] in dealing with their term J' (e) following the statement of their Lemma 5.1, by 
choosing p small enough. It is thus sufficient to prove that 



J{e) :=^E 



dudv [Xu+e - Xu ) (A„+e - Xy) g[ ) g 



tends to as e tends to 0. We now use the same method and notation as in Step 3 of the proof 
of Theorem 4.1 in [12]. In order to avoid repeating arguments from that proof, we only state and 
prove the new lemmas which are required. 

Step 1: translating Lemma 5.3 from JT^ . Using the fact that E < E [GJ\ < 6'^ (e), this lemma 
translates as: 

Lemma 14 Let k >2 be an integer. Then 

dudv < est ■ e6^ (e) . 



11^ 


|1 i\ 







This step and the next 4 steps are devoted to the Proof of lemma I14L We only need to show 
that for all i, j G {1,2}, 

\rij\'' dudv < est ■ e5^ {e) . (38) 



Recall the function K defined in [12 



K {u, v) := E [{Xu+e + Xu) {X,+e + X,)] 

= Q [u + e,v + e) + Q {u,v + e) + Q [u + e,v) + Q {u,v) . 

This is not to be confused with the usage of the letter K in previous sections, to which there will 
be made no reference in this proof; the same remark hold for the notation A borrowed again from 
[12] , and used below. 

To follow the proof in [12j, we need to prove the following items for some constants ci and C2'. 

1. ci5'^{u) < K{u,u) < C26^{u); 

2. K{u,v) < C25{u)6{v) ; 

3. A {u, v) := K (n, u) K {v, v) - K (u, vf > ciJ^ (^) p {v - u) . 

The lower bound in item 1 is our hypothesis (i). Also note that the hypothesis borrowed from 
Theorem [71 which says that 5^ (s,t) < i\t — implies, for s = 0, that 

6^{Q,u) = Qu<5'^{u). (39) 
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The upper bound in item 1 is a special case of item 2, which we now prove. We write, via Schwartz's 
inequahty and the fact that 5"^ is increasing, and thanks to p9p. 

K {u, v) <46{u + e)6{v + e). 

However, since (5^ is concave with 6 (0) = 0, we have 6^ {2u) /2u < 6^ (u) /u. Also, since we are in 
the set D^, u + e < 2u and v + e < 2v. Hence 

K{u,v) < 4d{2u)6{2v) 
< 86 [u) 6 {v) , 

which is item 2. By items 1 and 2, item 3 is a consequence of the existence of c > such that 

A {u,v) > cQuS'^ {u,v) 

= cQu (Qu + Qv-2Q (n, v)) . 

An elementary computation with A written using only Q shows that this last condition is equivalent 
to our hypothesis (ii). 

It will now be necessary to reestimate the components of the matrix A21 where we recall 



A2i[ll] 
A2i[12] 
A2i[21] 
A2i[22] 



E [(-'^»;+e + Xy) {Xu+s — X-u)] , 
E [{X-u+s + ^u) {Xv+s — X-v)] , 
E [(-'^^»;+e + Xy) — X^)] . 



Step 2: the term r\\. We have by the lower bound of item 1 above on K (n, n), 

1 



A2l[ll] 



c-st , . , 
< — |A2l[ll]|. 



To bound |A2i[ll]| above, we write 



|A2i[ll]| = |E + X„) - 

— Qu+e Qu 

< eQ (u) ju 

< e6^ (u) ju 

where we used the facts that is increasing and concave, and that < (5^ {u). Thus we have 

b{u) 



Fill < £ cst- 



u 



The result (j38|) for i = j = 1 now follows by the next lemma. 

Lemma 15 For every k > 2, there exists > such that for every e E (0, 1), 





6{u) 


k 


He) 




du < Ck£ 




u 




e 
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Proof of lemma 1151 Our hypothesis (iii) can be rewritten as 



6 (au) 



au 



< 



l + (a-l)6\ 6{u) 



a 

5{u) 



u 



u 



The concavity of 5 also imphes that 5 (u) /u is increasing. Thus we can write 



6{u) 



: OO 

j=0 



6{u) 



u 



du 



OO 



+1 



j=0 
£ {a 



5{e) 



He) 



k OO 



j=0 



The lemma will be proved if we can show that / (a) := \Ka^b\ a < 1 for some a > 1. We have 
/ (1) = and /'(I) = k{l — b) — 1. This last quantity is strictly positive for all A; > 2 as soon as 
h < 1/2. This finishes the proof of the lemma [T5l □ 



Step 3: the term ri2- We have 
ri2 = A21 [11] 



-K (n, v) 



a/ K {u, u) a (n, v) 



+ A21 [12] 



y^K iu,u) 



We saw in the previous step that IA21 [11]| = \Qu+e — Qu\ < est eS"^ (n) /u. For A21 [12], using the 
hypotheses on our increasing and concave functions, we calculate 

I A21 [12] I = |2 {Qu+e - Qu) +S'^{u + e,v + £)-6^ {u, v + e) + 6^ {u + e,v) - 5^ (n, v) \ 

< 2IA21 [11] I + £6^ {u + e,v + £)/{v-u) + eS^ {u + e,v) / {v - u - e) 

< 2 IA21 [11]| + £6^ {v -u) / (v -u) + £5'^ (v -u- e) / {v -u- e) 

< 2cst £6^ [u] /u + 2£6^ {v -u- £) / {v -u- £) . (40) 

The presence of the term —£ in the last expression above is slightly aggravating, and one would 
like to dispose of it. However, since {u,v) G D^, we have v — u > £'' ioi some p E (0, 1). Therefore 
V — u — £>£p — £> £^1"^ for £ small enough. Hence by using p/2 instead of p in the definition of 
in the current calculation, we can ignore the term —£ in the last displayed line above. Together 
with items 1, 2, and 3 above which enable us to control the terms K and A in ri2, we now have 



1^2! < est £ 
+ est £ 

= est £ 



S (u) 5 (v) 



+ 



5iu) 



5^ ju) 

u \6 (u) 6 (u) 5 (v — u) ' 6 (u) 5 {v — u) 
6^ (v-u) 6 (u) 
V — u 6 (u) 5 {v — u) 



6 ju) 6 jv) 6^ ju) 6 (v-u) 
uS {v — u) u5 {v — u) V — u 
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We may thus write 



ki2|'^ dudv < est 



5 (n) 6 {v) 



u5 {v — u) 



+ 



6^ (u) 



u6 {v — u) 



+ 



6 {v — u) 



V — u 



dudv. 



5{v—u) 



dudv is identical, after a trivial change of variables, to the one dealt 



u5(v—u) 



dudv is smaller than the first 



The last term //^ 

with in Step 2. Since 5 is increasing, second the term J jj^ 

k 

term JJ^ ^u5{v-u) dudv. Thus we only need to deal with that first term; it is more delicate than 
what we estimated in Step 2. 

We separate the integral over u at the intermediate point v/2. When u G [v/2,v — e], we use 
the estimate 

5{u) ^ 5 [v/2) ^Jiv)_ 
u ~ v/2 ~ V 

On the other hand when u G [e, w/2] we simply bound 1/5 [v — u) by 1/5 {v/2). Thus 

k 



< 



< 2^ 



5(u) 


5{v) 


u5 (v 


-u) 


dv 




2e 


J u=e 


dv 

2e 


5{v 


5{v/ 


r 


5{u) 


' u=e 


u 



dudv 



5 (n) 5 (v) 



u5 [v — u) 

k nv/2 



du + 

k 



dv 



5 (n) 5 {v) 



5{u) 



u 



du + 2 



u=v/2 
1 52 



du + 2 



S'' (e) J,.. 



f v' 


5{v) 






1 V=£ 


V 



V=£ 

2k 



u5 [v — u) 

k 

dv 



du 



V 



u=v/2 



1 



5 [v — u) 



du 



dv 



< est e 



5(e] 



here we used the concavity of 5 to imply that 5 (v) /5{v/2) < 2, and to obtain the last line, we 
used Lemma [12] for the first term in the previous line, and we used the fact that 5 is increasing and 
that V < 1, together again with Lemma [T5] for the second term in the previous line. This finishes 
the proof of (f38l) for ru. 



Step 4-' the term r2i. We have 



r2i = A21 [21] 



\/K {u,u) 



and similarly to the previous step. 



IA21 [21]| = \Q {u + e,v + e) - Q {u + £,v) + Q {u,v + e) - Q {u,v)\ 

= |2 (Q^+e - Qv) +5'^{u + e, v)-5'^{u + e,v + e) + 5^ (n, v) - 5^ (n, v + e)\ 

, 5'^ (u + e,v) 5"^ (n, v) 

< 2 A21 11 — +£ ^ ' ' 

V — u — e V — u 

< 2cst £(5^ (u) /u + AeS"^ (v - u) / {v - u) , 
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which is the same expression as in ()40p . Hence with the lower bound of Item 1 on X (u, u) we have 



\r2i\ dudv < est e 



est e" 



-De 



5(u) 



6{u) 



+ 



+ 



6'^ {v — u) 



{v — u) 6 (n) 



6^ (u) 



u5 [v — u) 



dudv 



dudv. 



This is bounded above by the expression obtained as an upper bound in Step 3 for JJp ^12!^ dudv, 
which finishes the proof of ()38p for r2i. 



Step 5: the term r22- Here we have 



r22 = A21 [21] 



-K (n, v) 



\J K {u, u) A (n, v) 
We have ah'eady seen in the previous step that 



+ A21 [22] 



VK {u,u) 
\/A(n,?;)' 



IA21 [21] I < est e 
Moreover, we have, as in Step 2, 



5"^ {u) ^5'^{v- u) 



v — u 



IA21 [22] I = \Q^+s-Qv\ <cste 



52 (v) 



Thus using the bounds in items 1, 2, and 3, 



|?^22| < est e 
= est e 



(52 (n) 5'^{v-u)\ 5{u)5{v) 5^ {v) 5{u) 



+ 



u 



+ 



V — u J 6"^ (u) 6 {v — u) v 6 (u) 5 {v — u) 



6{u)6{v) 6{v)6{v-u) 6'^{v) 



+ 



u6 [v — u) 5 (n) {v — u) v6 {v — u) 



Of the last three terms, the first term was already treated in Step 3, the second is, up to a change 
of variable, identical to the first, and the third is smaller than f, ^ which was also treated in 

' ' ud(v—u) 

Step 3. Thus (j38p is proved for r22, which finishes the entire proof of Lemma □ 
Step 6: translating Lemma 5.4 from 112]. We will prove the following result 
Lemma 16 For all j e {0, 1, • • • , (m - 1) /2}, 

jj |E [Z3Z4]|'"-^^' dudv < est e(52(™-2i) (e) . 

Proof of Lemma 1161 As in [12] , we have 

|E [Z^Z^ir-^^ < est |E [GsG^ir-^^ + est |E [TaTJI'^-^J' . 

The required estimate for the term corresponding to |E [r3r4]|™'~^-' follows by Schwartz's inequality 
and Lemma [HI For the term corresponding to |E [G3G4]|'"~^"' , we recognize that E [G3G4] is the 
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negative planar increment 0^(n, defined in (jlOp . Thus the corresponding term was already 
considered in the proof of Theorem d?]). More specifically, up to the factor e'^5~'^^ (e), we now have 
to estimated the same integral as in Step 2 of that theorem's proof: see expression (|22p for the 
term we called Jj,OD- This means that 

I IE [GsCr-^^ dudv < ^^J,,OD 

<e2|/i|(OL>)<52(™-2i-i) (e). 

Our hypotheses borrowed from Theorem ([7|) that |^| [OD) < est ei/™"^ and that 5^ (e) = o (r^/*^^™")) 
now imply that the above is <C eJ^^™-"^-?) (e), concluding the lemma's proof. □ 

Step 7. Conelusion. The remainder of the proof of the theorem is to check that Lemmas O and [16] 
do imply the claim of the theorem; this is done exactly as in Steps 3 and 4 of the proof of Theorem 
4.1 in [12]. Since such a task is only bookkeeping, we omit it, concluding the proof of Theorem 1121 



References 

[1] Alos, E., Mazet, O., Nualart, D.: Stochastie ealeulus with respect to Gaussian processes. Ann. 
Probab. 29, 766-801 (1999) 

[2] Bertoin, J.: Sur une integrale pour les processus a a-variation bornee. Ann. Probab., 
17(4):1521-1535, 1989. 

[3] Burdzy, K., Swanson, J.: A change of variable formula with ltd correction term. Preprint, 
available at http:/ /arxiv.org/PS_cache/arxiv/pdf/0802/0802. 3356v3.pdf. 

[4] Bruneau, M.: Variation totale d'une fonction. Lecture Notes in Mathematics, Vol. 413. 
Springer- Verlag, Berlin-New York, 1974. 

[5] Dudley, R. M., Norvaisa, R.: Differentiability of six operators on nonsmooth functions and 
$p$-variation. With the collaboration of Jinghua Qian. Lecture Notes in Mathematics, 1703. 
Springer- Verlag, Berlin, 1999. 

[6] Errami, M., Russo, F.: Covariation de convolution de martingales. C. R. Acad. Sci. Paris 
Ser. I Math. 326 (1998), no. 5, 601-606. 

[7] Flandoli, F., Russo, F.: Generalized stochastic integration and stochastic ODE's. Annals of 
Probability. Vol. 30, No 1, 270-292 (2002). 

[8] Follmer, H.: Calcul d'lto sans probabilites. In Seminaire de Probabilites, XV (Univ. Stras- 
bourg, Strasbourg, 1979/1980) (French), volume 850 of Lecture Notes in Math., pages 143-150. 
Springer, Berlin, 1981. 

[9] Gradinaru, M., Russo F., Vallois, P.: Generalized covariations, local time and Stratonovich 
Ito's formula for fractional Brownian motion with Hurst index H > j. Ann. Pobab., 31 (4), 
1772-1820, 2003. 



38 



[10] Errami, M., Russo, F.: n-covariation, generalized Dirichlet processes and calculus with respect 
to finite cubic variation processes. Stochastic Processes Appl. , 104, 259-299 (2003). 

[11] Gradinaru, M., Nourdin, I.: Approximation at first and second order of m- order integrals of 
the fractional Brownian motion and of certain semimartingales. Electron. J. Probab. 8, no. 
18, 26 pp (2003). 

[12] Gradinaru, M.; Nourdin, I.; Russo, F.; Vallois, P.: m-order integrals and Ito's formula for 
non-semimartingale processes; the case of a fractional Brownian motion with any Hurst index. 
Ann. Inst. H. Poincare Probab. Statist. 41, 781-806 (2005). 

[13] Hu, Y., Nualart,D., Song, J.: Fractional martingales and characterization of the fractional 
Brownian motion. Preprint, http://arxiv.org/abs/0711.1313. 

[14] Lyons, T., Qian, Z.: System control and rough paths. Oxford Mathematical Monographs. 
Oxford University Press, Oxford, 2002. 

[15] Mocioalca, O., Viens, F.: Skorohod integration and stochastic calculus beyond the fractional 
Brownian scale. Journal of Functional Analysis, 222 no. 2, 385-434 (2004). 

[16] Nourdin, I.: A change of variable formula for the 2D fractional Brownian motion of Hurst 
index bigger or equal to 1/4. J. Funct. Anal. 256, 2304-2320 (2009). 

[17] Nourdin, I., Nualart D., Tudor, C.. Central and non-central limit theorems for 
weighted power variations of fractional Brownian motion. Preprint, available at 
http://arxiv.org/PS-cache/arxiv/pdf/0710/0710.5639vl.pdfi 

[18] Nualart, D., Ortiz-Latorre, S.: Central limit theorems for multiple stochastic integrals and 
Malliavin calculus. Stochastic Process. Appl., 118(4), 614-628 (2008). 

[19] Nualart, D.: The Malliavin calculus and related topics. Second edition. Probability and its 
Applications. Springer- Verlag, 2006. 

[20] Rogers, L. C. G., Walsh, J. B.: The exact A/S-variation of a process arising from Brownian 
motion. Stochastics Stochastics Rep. 51 (1994), no. 3-4, 267-291. 

[21] Russo, F., Vallois, P.: Stochastic calculus with respect to a finite variation process. Stochastics 
and Stochastic Reports 70, 1-40 (2000). 

[22] Russo, F., Vallois, P.: The generalized covariation process and ltd formula. Stochastic Process. 
Appl. 59, no. 1, 81-104 (1995). 

[23] Russo, F., Vallois, P.: Elements of stochastic calculus via regularizations. Seminaire de Prob- 
abilites XL, Lect. Notes Math. 1899, Berhn Heidelberg New- York, Springer, 147-186 (2007). 

[24] Viens, F., Vizcarra, A.. Supremum Concentration Inequality and Modulus of Continuity for 
Sub-nth Chaos Processes. J. Funct. Analysis 248, 1-26 (2007). 



39 



